Fractional Talbot effect in phase space: 
A compact summation formula 
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Abstract: A phase space description of the fractional Talbot effect, 
occurring in a one-dimensional Fresnel diffraction from a periodic grat- 
ing, is presented. Using the phase space formalism a compact summa- 
tion formula for the Wigner function at rational multiples of the Talbot 
distance is derived. The summation formula shows that the fractional 
Talbot image in the phase space is generated by a finite sum of spatially 
displaced Wigner functions of the source field. 
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It is the purpose of this communication to investigate a one dimensional Fresnel diffrac- 
tion from a periodic grating and the corresponding fractional Talbot effect using the 
Wigner phase space distribution function. 



We use the following definition of the Wigner distribution function |j for the 
held amplitude E(x) along the axis x perpendicular to the propagation direction: 



W(x,u) = ±- JdyE*(x + y/2)e m yE(x-y/2), (1) 

where u has the meaning of a spatial frequency. We consider a monochromatic plane 
wave characterized by its electric held £ (x; z) = e tkz E(x; z) propagated paraxially along 
the z-axis. At z = 0, where the field starts to propagate, there is an infinite one- 
dimensional periodic grating with transmittance t(x) = t(x + a). The electric held 
amplitude after the passage through the grating can be expanded into the Fourier series: 

oo 

E(x;0)=t(x)E o = tne 27rmx/a , (2) 



n— — oo 



where a is the period of the grating and Eq is a constant amplitude of the incident wave. 
The Wigner distribution function of this source held is given by: 



W(x,u;0) = J2 K\ 2 S(u-nu ) 



+ tntn' exp[2ni(n — n')x / a\5(u — — uq). (3) 

Two different types of contributions can be distinguished in the above formula. The 
first one is a set of parallel positive density stripes separated by the spatial frequency 
uq = 2ir/a, given by the reciprocal grid spacing. Each of them is generated by a sepa- 
rate Fourier component of the source field. Coherence between the Fourier components 
results in nonpositive oscillatory terms of the Wigner function, located precisely in the 
middle between the contributing frequencies. The nonpositive interference terms arc 
a consequence of the linear superposition principle and the bilinear character of the 
Wigner function. 

The paraxial propagation of the field through a distance z is described in the 
phase space as the following simple transformation of the Wigner distribution function 
1: 

W(x,u;z) = W(x- -^u,u;0). (4) 
2n 

This formula applied to the Wigner function of the periodic field gives the following 
result: 



W(x,u;z) = \t n \ 2 S(u — nuo) 



n— — oo 



+ ^2 tnt *n' exp[27ri(n - n')x/a - 2ixi(B n - 6 n >)]8(u — — u ). (5) 

In the course of propagation the interference terms acquire additional phase shifts given 

by 



z 



,2 



= — nr, (6) 
z T 



where z? = a 2 /2A is the Talbot distance. It is straightforward to see that for inte- 
ger multiples of the Talbot distance the original Wigner function of the input field is 
reproduced. 



At intermediate distances, the phase shifts 9 n play a nontrivial role and the 
structure of the observed Fresnel images becomes more complex. Nevertheless, they 
exhibit an interesting regular behavior at rational multiples of the Talbot distance. We 
will now discuss this effect in terms of the Wigner distribution function. Let us denote 
z/zt = p/q, where p and q are coprime integers. The main complication in Eq. (|^) are 
phase factors 9 n which depend quadratically on n. We will simplify it with the help of 
an observation used in the studies of quantum wave packets dynamics Q : the exponent 
exp(— 2ni9 n ) is periodic in n with the period I = q/4 if q is a multiple of 4 and I = q 
otherwise. The quadratic phase factor can therefore be represented as a finite Fourier 
sum: 



exp(— 2iri0 n ) 



l-i 

E 

s=0 



a s exp(—2irisn/l) 



(7) 



with certain coefficients a s , which have been analyzed in detail in Ref. || . Substitution 
of the above expression yields: 



W(x,u;pz T /q) 



1-1 

s,s'=0 n,n' — — oo 

2-nin' 



E27rm / sa\ 
t n t n , exp — — [ x ~ —) 
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\ U Uq 



(8) 



The interference terms of the Wigner function can now be interpreted as generated 
by pairs of Fourier components of the source field shifted in the position space by 
rational fractions sa/l, where a = 0,1, ... ,1 — 1. A simple rearrangement of the exponent 
arguments allows one to represent the sum over n and n' solely in terms of the source 
field Wigner function: 



l-i 



W(x,u;pzt /q) = a s a*j exp[iu(s' — s)a/l]W I x — 



s,s'=0 



Q + s')a 
21 : 



u-0 



(9) 



This is a compact summation formula for the Wigner function at rational multiples of 
the Talbot distance. It shows, quite surprisingly, that W(x, u;pzT/q) is simply given by 
a finite sum of spatially displaced Wigner functions of the source field, with some phase 
factors. 

Integration of the derived expression for the Wigner function over u yields the 
known formula for the field intensity distribution in the observation plane: 



\E(x;pz T /q)\' 



du W {x , u; pzt / q) 



i-i 



a s t{x — sa/l; 0) 



(10) 



which shows that the observed Fresnel image can be represented as generated by a finite 
sum of shifted source field amplitudes 0. A detailed discussion of fractional Talbot 
images can be found in Ref. 0. 

We have discussed a one dimensional Fresnel diffraction from a periodic grating 
using the phase space formalism. We have seen in this picture, that the whole propaga- 
tion of the field is encoded in the interference terms of the Wigner function, generated 
by coherence between Fourier components of the source wave. The variety of the Fresnel 
images is simply a result of a complicated interplay between phases of these interference 



terms. The summation formula derived in this paper demonstrates that the fractional 
Talbot effect can be also understood as an interference between components shifted in 
the position space by fractions of the pattern period. 

Let us finally note a close analogy between the Talbot effect and the dynam- 
ics of various nonlinear quantum systems. Dynamical phase factors with a quadratic 
dependence on the summation index are a common feature of many problems. These 
phase factors can be easily revealed in the equation of motion for the quantum Wigner 
function in the (x,p) phase space: 



where ip n are the amplitudes of the initial wave function ?p(x; 0) projected onto the 
energy eigenstate tp n (x), and where 



are the cross Wigner functions of the eigenstates. For many systems, their eigenergies E n 
are given by a quadratic polynomial of n. This can be either an approximate dependence, 
as it is for the Rydberg electron wave packets M, or a strict one, as it happens for an 
electromagnetic field mode in the Kerr medium [[7[ and a particle in the infinite square 
well. In all these cases, the state of the system turns out to be highly regular for times 
equal to rational multiples of the characteristic revival time, at which the system is 
again in its initial state. This regularity consists in a formation of a finite superposition 
of copies of the initial state generated by a simple transformation. For the wave packets 
dynamics, the components are spatially shifted ||; for the Kerr medium, the components 
are rotated in the phase space around its origin j^] . 

In particular, a close similarity exists between the Talbot effect and the frac- 
tional revivals in the infinite square well. The dynamics of the latter system is equivalent 
to an evolution in a free space of an infinite periodic sequence of displaced initial wave 
functions Furthermore, the propagation of electromagnetic fields in the paraxial 

approximation is governed by the same wave equation as the quantum time evolution. 
Therefore, Fresnel diffraction of classical electromagnetic waves on periodic gratings is 
a strict counterpart of the quantum dynamics of a confined particle. 

Acknowledgments 

The authors have benefited from discussions with J. H. Eberly and from a stimulating 
talk given in Ulm by B. Rohwedder. This work has been partially supported by the 
Polish KBN grant 2 P03B 006 11 and the Deutsche Forschungsgemeinshaft. 




(11) 




(12) 



